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available on astroph given the rapid progress in the field of GRBs. You are free to refer 
to this paper in your own paper. However, we do place restrictions on any dissemina- 
tion in the popular media. The article is under embargo until it is published. For fur- 
ther enquiries, please contact Shri Kulkarni (srk@astro. caltech. edu) or Joshua Bloom 
(j sbOastro . caltech . edu). 

Cosmic gamma-ray bursts (GRBs) have been firmly established as one of the 
most powerful phenomena in the Universe, releasing electromagnetic energy 
approaching the rest-mass energy of a neutron star in a few seconds. The 
two currently popular models for GRB progenitors are the coalescence of two 
compact objects (such as neutron stars or black holes) or collapse of a massive 
star. An unavoidable consequence of the latter model is that a bright super- 
novae should accompany the GRB. The emission from this supernova competes 
with the much brighter afterglow produced by the relativistic shock that gives 
rise to the GRB itself. Here we present evidence for an unusual light curve 
for GRB 980326 based on new optical observations. The transient brightened 
~ 3 weeks after the burst to a flux sixty times larger than that extrapolated 
from the rapid decay seen at early time. Furthermore, the spectrum changed 
dramatically and became extremely red. We argue that the new source is the 
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underlying supernova. If our hypothesis is true then this would be the first 
evidence for a supernova connection with GRBs at cosmological distances. We 
suggest that GRBs with long durations are associated with death of massive 
stars. 

The origin of GRBs remained elusive for a period of nearly three decades after their 
discovery 1 . Beginning in 1997, however, the prompt localization of GRBs by the Italian- 
Dutch satellite BeppoSAX 2 and the All Sky Monitor 3 on board the X-ray Timing Explorer 
led to the discovery of the GRB afterglow phenomenon - emission at lower energies: X- 
ray 4 , optical 5 , and radio 6 . 

The persistence of the afterglow emission (days at X-ray wavelengths, weeks to months 
at optical wavelengths, months to a year at radio wavelengths) enabled astronomers to 
carry out detailed observations which led to fundamental advances in our understanding 
of these sources: (1) the demonstration that GRBs are at cosmological distances 7 ; (2) the 
proof that these sources expand with relativistic speeds 6 ; and (3) the realization that the 
electromagnetic energy released in these objects exceeds that in supernovae 8 and, in some 
cases, the released energy is comparable to the rest mass energy of a neutron star 9 ' 10 ' 11 ' 12 . 

Despite these advances, we are still largely in the dark about the nature of the GRB 
progenitors. Though there are a number of models for their origin, the currently popular 
models involve the formation of black holes resulting from either the coalescence of neutron 
stars 13 ' 14 ' 15 or the death of massive stars 16 ' 17 . The small offsets of GRBs with respect to 
their host galaxies and the association of GRBs with dusty regions and star-formation 
regions favors the latter, the so-called hypernova scenario 17 . However, this evidence is 
indirect and also limited by the small number of well-studied GRBs. 

The most direct evidence for a massive star origin would be the observation of a 
supernova coincident with a GRB. Here we present observations of GRB 980326 and argue 
for the presence of such an underlying supernova. If our conclusions are correct then 
the implication is that at least some fraction of GRBs, perhaps the entire class of long 
duration GRBs, represent the end point of the most massive stars. Furthermore, if the 
association 18 ' 19 of GRB 980425 with a bright supernova in a nearby galaxy holds, then the 
apparent 7-ray luminosity of GRBs ranges over six orders of magnitude. 

The unusual optical afterglow 

Following the localization of GRB 980326 by BeppoSAX (ref. 20), Groot et al. 21 
quickly identified the optical afterglow. Our optical follow-up program began at the Keck 
Observatory, just 10 hr after the burst. A log of these observations is given in Table 1. 

In Figure 2 we present our .R-band photometry along with those reported in the lit- 
erature. Restricting to data taken within the first month of the burst reported in the 
literature 21 ' 22 and via the "GRB Coordinates Network" (GCN) (ref. 23), we find a charac- 
teristic power law decay in the flux versus time followed followed by an apparent flattening. 
The usual interpretation is that the decaying flux is the afterglow emission, while the con- 
stant flux reflects the presence of the host galaxy. Indeed, earlier 24 we attributed the entire 
observed flux on April 17th to the host galaxy. 
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However, to our surprise, our more recent observations (first carried out nine months 
after the GRB event) showed no galaxy at the position of the optical transient (OT); see 
Figure 1. We estimate a 2-a upper limit of R > 27.3 magnitude (see Table 1). This 
is almost a factor of 10 less flux than that reported from our April 17th detection. A 
secure conclusion is that the presumed host galaxy of GRB 980326, assuming the GRB 
was coincident with the host (as appears to be the case for all other well-studied GRBs to 
date) is fainter than R ~ 27 magnitude. This conclusion is not alarming since such faint 
(or fainter) galaxies are indeed expected from studies 25 ' 26 of the properties of cosmological 
GRB host galaxies. 

Having established that the host galaxy of GRB 980326 is faint, we are forced to 
conclude that the OT did not continue the rapid decay it exhibited initially. Instead, we 
find two phases of the light curve (Figure 2): a steeply declining initial phase (t ^ 5 d) and 
a subsequent rebrightening phase (t ~ 3-4 weeks). Following the rebrightening, the source 
appears to have faded away to an undetectable level by the time of our next observation 
(9 months after the burst). 

In previously studied bursts the optical afterglow emission has been modeled by a 
power law function, flux oc t a ; here, t is time since the burst and a the power law index. 
In some bursts, at early times (t less than a day or so), significant deviations have been 
seen, e.g., GRB 970508 (ref. 27). At late times, in some bursts, deviations manifest as 
steepening (i.e. a becoming smaller) of light curves, e.g., GRB 990123 (ref. 11). 

It is against this backdrop of the observed afterglow phenomenology that we now 
analyze the light curve in Figure 2. The declining phase cannot be fit by a simple power 
law (x 2 = 72 for 9 degrees of freedom). From Figure 2 it is clear that the flux already starts 
flattening by day 3. Restricting the analysis to the first two days, we obtain a = —2.0 ±0.1, 
consistent with previous analysis 21 . 

Such power law decays are usually interpreted as arising from electrons shocked by the 
explosive debris sweeping up the ambient medium 28 ' 29 ' 30 ' 31 . Assuming that the electrons 
behind the shock are accelerated to a power-law differential energy distribution with index 
— p, on general grounds 32 we expect that the afterglow flux, f u (t) oc t a v^\ here f u (t) is 
the flux at frequency v and time t. The value of a and /3 depend on p, the geometry 
of the emitting surface 33 (spherical versus collimation) and the radial distribution of the 
circumburst medium 34 . 

From our spectroscopic observations of March 29th (Figure 3) we find (3 = —0.8 ±0.4. 
This combination of (a, j3) is similar to the (a = —2.05 ± 0.04, (3 = —1.20 ± 0.25) seen 
in GRB 980519 (ref. 35) and can be reasonably interpreted 36 as arising from a standard 
p ~ 2.2 shock with a jet-like emitting surface. Alternatively, the emission could arise in 
ap~3 shock propagating in a circumburst medium 34 whose density falls as the inverse 
square of the distance from the explosion site. 

A new transient source 

We now discuss the bright source seen in the rebrightening phase (corresponding to 
observations of April 17 and April 23). This source is ~ 60 times brighter than that 
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extrapolated from the rapidly declining afterglow. The magnitude of this excess and the 
late timescale of rebrightening has never been reported before. 

We first establish the reality of the source. As noted in the legend to Figure 1 the 
source is consistently detected in three separate images of April 17. In the summed image, 
the source is detected at 4.6-cr (chance probability of 2 x 10~ 6 ). We note that all other ob- 
jects in the field at this flux level are reliably detected in our deeper December 18th image. 
Next, the source is clearly detected in the spectrum obtained on April 23rd (Figure 3) and 
at the same position as that of the OT. Finally, we note that the source in the April 17th 
image is coincident with the position of the OT in the image of March 27th to within the 
expected astrometric error, 0.04 ± 0.18 arcsecond. 

We conclude that there was indeed a source at the position of the OT which bright- 
ened three weeks after the burst and subsequently faded to undetectable levels. We now 
investigate possible explanations for this source. 

The simplest picture is that the rebrightening phase is due to a rebrightening of 
the optical afterglow itself. As noted earlier, this would be unprecedented in both the 
magnitude of the rebrightening and the epoch of rebrightening. Piro et al. have recently 
suggested that the doubling of the X-ray flux of GRB 970508 three days after the GRB 
event arises from the relativistic shell running into a dense gas cloud. Such an explanation 
for the GRB 980326 light curve would require a large dense region, with a size comparable 
to the timescale of rebrightening, about T x 10 light days (~ 0.1 pc) and located at a 
distance Y 2 c x 20 days (~ 1 pc) from the explosion site. Here, Y is the bulk Lorentz factor 
of the shock and is expected to be order unity three weeks after the burst. Panaitescu 
et al. 38 suggest the rebrightening of GRB 970508 may be due to a shock refreshment - 
delayed energy injection by the extremely long-lived central engine that produced the GRB. 
In both these models, the expected spectrum would be the typical synchrotron spectrum, 
flux f v oc v~ x (or flatter). The very red spectrum of April 21st (Figure 3) allows us to 
essentially rule out a synchrotron origin for the rebrightening phase of GRB 980326. 

Alternatively, the GRB could have occurred in a dusty region and the afterglow would 
rebrighten 39 as the dust is sublimated by the afterglow. However, the observed spectral 
evolution from a relatively blue spectrum (March 29) to red (April 23) moves in a direction 
opposite to that expected in this scenario. 

The supernova interpretation 

We advance the hypothesis that the new source is due to an underlying supernova 
(SN) revealed only after the afterglow emission has vanished. Woosley and collaborators 
(see ref. 40,16 and references therein) have pioneered the "collapsar" model in which GRBs 
arise from the death of massive stars - stars which produce black hole remnants rather than 
neutron stars. In this model, the iron core of a massive star collapses to a black hole and 
releases up to a few xlO 52 erg of kinetic energy. Some fraction of this energy is expected 
to emerge in the form of a jet with little entrained matter; bursts of gamma-rays result 
from internal shocks in this jet. The remaining energy is absorbed by the star, causing it 
to explode and thereby produce an energetic supernova. 
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Thus in this model, the total light curve has two distinct contributions: a power-law 
decaying afterglow component and emission from the underlying supernova. In Figure 2 
we present the light curve expected in this model and use the light curve of the well 
observed 18 ' 41 SN 1998bw as a template for the supernova contribution. We find the R- 
band and /-band data consistent with a bright supernova at z ~ 1. 

The very red spectrum of the source on April 23 finds a natural explanation in the su- 
pernova hypothesis. From theoretical 40 and phenomenological 19 grounds we expect GRBs 
to arise from massive stars which have lost their hydrogen envelope, i.e. Type Ibc super- 
novae. At low redshifts, all Type I supernovae are observed to exhibit a strong UV deficit 
relative to the blackbody fit to their spectra. This deficit is due to absorption by prominent 
atomic resonance lines starting below ~ 3900 A. Below A c ~ 2900 A we expect to see very 
little flux. In the near UV range (3000-4000 A) type I SNe spectra have a red appearance. 
Approximating the flux by a power law (/„ oc the power law index (depending on 
the wavelength range chosen) is —3 or even smaller; see ref. 42 for a UV spectrum of a 
type la SN. Fitting the spectrum of Figure 3 to a power law we obtain f3 = — 2.8 ± 0.3. In 
this interpretation the redshift of the source is z ~ 1. A smaller redshift would lead to a 
larger (3. A larger redshift would substantially suppress the light in the observed R band 
(which covers the wavelength range 5800-7380 A). Indeed, from Figure 1, we can deduce 
that z < 1.6. 

We do not know a priori the spectrum and light curve of a supernova accompanying 
GRBs. However, we have used the light curve of SN 1998bw because it is a very well studied 
Type Ibc SN with a possible association with GRB 980425. If SN 1998bw is indeed an 
appropriate template for a supernova associated with a GRB then, as Figure 2 illustrates, 
the likely redshift of GRB 980326 is z ~ 1 and, as discussed above, the red spectrum of 
April 23 similarly suggests z ~ 1. Given the low signal-to- noise ratio of the April 23rd 
spectrum and the expected line broadening due to high photospheric velocity we do not, 
as seems to be the case, expect to see any spectral features. 

Independently, from the absence of strong spectral breaks in our spectrum of the OT 
we can firmly place zot <, 2.3. This constraint is consistent with our deduction that 
z <^ 1.6 (see above). Thus from a variety of accounts we find a plausible redshift of around 
unity for GRB 980326. Such a redshift is not entirely unexpected. Indeed, we note that five 
out of eight spectroscopically confirmed redshifts of GRBs lie in the range 0.7 < z < 1.1. 

Implications of the supernova connection 

The GRBs localized by BeppoSAX belong to a class of long duration GRBs. The jet 
in a collapsar model takes many seconds to penetrate the star, and therefore the collapsar 
model is unlikely to account 43 for the class of short duration (less than a few seconds) 
GRBs. If we accept the SN interpretation for GRB 980326, a long duration (5 s) GRB 
then it is only reasonable to posit that all other long duration GRBs are also associated 
with SNe. In what way can this assertion be tested observationally? 

The evidence for an underlying SN can come in two ways. First, is the direct evidence 
for an accompanying SN seen in the light curve at timescales comparable to the time for 
SNe to peak, ~ 20(1 + z) days. However, in our opinion three conditions must be satisfied 
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in order to see the underlying SN even when one was present. (1) The GRB afterglow 
should decline rapidly, otherwise the SN will remain overpowered by the afterglow for all 
epochs. (2) Given the strong UV absorption (discussed above), only GRBs with redshift 
z ^ 1.6 have an observable SN component in the optical band. (3) The host must be 
dimmer than the peak magnitude of the SN (My ~ —19.5). The last requirement is not 
needed if the GRB can be resolved from the host (e.g. with HST). Finally one caveat 
is worth noting: the peak magnitudes of Type Ibc SNe are not constant (unlike those of 
Type la) and can vary from —16 mag to a maximum of —19.5 mag; see ref. 44. We have 
investigated the small sample of GRBs with adequate long-term follow up and conclude 
that perhaps only GRB 980519 satisfies the first and the third observational conditions for 
SN detection; the redshift of this GRB is unfortunately unknown. 

The second method is an indirect method to tie GRBs and SNe. The dynamics of 
the relativistic blast wave is strongly affected by the distribution of circumstellar matter. 
Chevalier and Li 34 note that massive stars, through their active winds, leave a circumstellar 
medium with density falling as the inverse square of the distance from the star. One 
expects smaller a for GRBs exploding such a circumstellar medium. In this framework, 
GRB afterglows which decline rapidly and are at modest redshifts will again be prime 
targets to search for the underlying SN. 

In conclusion we note that it is not possible to firmly demonstrate on purely obser- 
vational grounds that all long duration GRBs can be explained by the collapsar model. 
However, we strongly urge sensitive observations especially at longer wavelengths (to avoid 
the UV cutoff of SNe) for GRBs satisfying the above three conditions. If the proposed 
hypothesis is correct then the light curves and the spectra would show the behavior shown 
and discussed in Figure 2 and Figure 3. 

We end with a discussion of one interesting point. The total energy release in 7-rays 
of GRB 980326 was E 7 = 3x 10 51 /j erg where fj is the fractional solid angle of the jet 
(if any); here we have used the measured fluence 21 and assumed z ~ 1 (H = 65 km s _1 
Mpc" 2 , O = 0.2, A = 0). If this GRB was beamed then E 1 ~ 10 49 erg. Curiously 
enough, this rather small energy requirement places GRB 980326 as close in energetics to 
GRB 980425 {E 1 = 6 x 10 46 erg, ref. 18) as to the classic gamma-ray bursts {E 1 ^6x 10 51 
erg). 
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Table 1. Keck II Optical Observations^ of GRB 980326* 



Date" 


Band/ 


Int. Time 


Seeing 


Magnitude^ 


Observers 


(UT) 


Grating (sec) 


(FWHM) 






Mar 97 1^ 


R 


240 


n" 74 


oi ok -i- n 

zi.zu in vj.vjo 


A VP DPT AflR 


Mar 98 95 


R 


240 


u .uu 




WC Af) TiQ QAQ 

no, -rvu, uo, oiT-O 


Mar 2Q 97 


300 


3600 




94 45 + n s c 


HS AD DS SAS 


Mar 30.24 


R 


900 


0".93 


24.80 ±0.15 


SP, BG, RK, IH 


Apr 17.25 


R 


900 


0"82 


25.34 ± 0.33 d 


PC, JB 


Apr 23.83 


300 


5400 




24.9±0.3 C 


SGD, SCO 


Dec 18.50 


R 


2400 


0".74 


> 27.3 e 


SRK, JSB, MvK 


Dec 18.54 


I 


2100 


0".74 


> 25.3 e 


SRK, JSB, MvK 


Mar 24 


I 


5450 


0".80 


> 26.6 e 


SRK, JSB 



Notes: 

f We used the Keck II 10-m Telescope 2,048 x 2,048 pixel CCD (charged coupled device) 
Low- Resolution Imaging Spectrometer 45 (LRIS) for imaging and spectroscopy of the GRB 
field. 

* The epoch of GRB 980326 is March 26.888, 1998 (ref. 20). 

(a) Mean epoch of the image. The year is 1998 for all images except for that on March 24 
for which it is 1999. 

(b) Photometric Calibration. The absolute zero-point of the R (effective wavelength 46 
of A e ff ~ 6588 A) and /-bands (A e ff ~ 8060 A) were calibrated to the standard Cousins 
bandpass using standard-stars in the field SA98 (ref. 47) and assuming the standard atmo- 
spheric correction on Mauna Kea (0.1 mag and 0.06 mag per unit airmass, respectively). 
The estimated statistical error on the absolute zero-point is 0.01 mag. We estimate the 
systematic error (due to lack of inclusion of color term) to be less than 0.1 mag. We prop- 
agated all photometry to the absolute zero-point derived in the first epoch of observation 
using 8 "secondary" stars which were detected with high signal-to-noise ratio, unsaturated, 
near to the transient, and common to every epoch; the typical uncertainty in the zero-point 
propagation is 0.01 mag. Thus any systematic error in our absolute zero-point will not 
affect the conclusions based on relative flux. The uncertainties quoted in the Table contain 
all known sources of error (aperture correction, etc.) The calibrated magnitudes of the 
secondary stars reported in Groot et al. 21 agree to within the measurement errors. 

(c) Spectrophotometric measurement. The flux in /^Jy is determined at 6588 A, the central 
wavelength of the R c band; the conversion to magnitude assumes mag equal to 3020 Jy 
(ref. 46) . The spectrophotometric magnitudes are relative to a bright star that was on the 
slit (for which we have obtained independent photometry from our images). 

(d) Photometry of the faint source. Since the transient was not detected to significantly 
fainter levels in later epochs it is safe to assume that the April 17 detection was that of a 
point-source (and not an extended galaxy as we had earlier 24 believed). To maximize the 
signal-to-noise ratio we choose to measure the photometry in an aperture radius equal to 
the FWHM of the seeing and correct for the missing flux outside the aperture by using 
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the radial flux profiles of bright isolated stars in the image. The determination of the 
optimum sky level (from which we subtract the total flux in the aperture) is not well- 
defined. We estimate the systematic uncertainty introduced by the uncertainty in the sky 
level as 0.25 mag. The statistical uncertainty (weighted mean over different background 
determinations) of the flux was 0.22 mag. Thus we quote the quadrature sum of the 
statistical and systematic uncertainty of 0.33 mag. 

(e) Upper-limits. On 1998 Dec 18 UT and 1999 March 24 UT there was no detectable flux 
above the background at the position of the optical transient. We centered 1000 apertures 
randomly in our image (approximately 1800 x 2048 pixels in size) and performed weighted 
aperture photometry with a local determination of sky background and recorded the counts 
("DN") above background at each location. The flux contribution from an individual 
pixel, some radius r from the center of the aperture, to the total flux was weighted by a 
Gaussian with a radial width FWHM equal to the seeing. A histogram of the resulting 
flux was constructed. This histogram was decomposed into two components — a Gaussian 
with median near zero DN and a long tail of positive DN corresponding to actual source 
detections. We fit a Gaussian to the zero-median component, iteratively rejecting outlier 
aperture fluxes. Based on the photometric zero-point and using isolated point sources in 
the image for aperture corrections, we computed the relationship between DN within the 
weighted aperture and the total magnitude. In the Table we quote an upper limit (95%- 
confidence level corresponding to 2-cr of the Gaussian fit) at the position of the optical 
transient. 
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Figure 1. 

Images of the field of GRB 980326 at three epochs. Each images shows a 54" x 
54" region centered on the optical transient (labeled "OT"). In all the images, the local 
background has been subtracted by a median filter and the resulting image smoothed (with 
a two-dimensional Gaussian with a = 0".23). An unrelated faint source "f" in the field is 
noted for comparison of the relative limiting flux between the three epochs: it is marginally 
detected (at the ~ 2-cr level) on March 27 and April 17 but well detected on December 18. 
In contrast the OT is brighter and better detected (at the 4.6-a level, see text) on April 
17 but clearly not detected to fainter levels on December 18 (R > 27.3; see Table 1). 
Methodology. In keeping with standard practice, our April 17th observations consisted 
of three separate 300-s observations (dithered by 5 arcseconds). Visual inspection of the 
three frames reveals a faint source near the position of the optical transient. In no frames 
was did a diffraction spike of the nearby bright star "A" overlap the OT position. Also, 
there were no apparent cosmic-ray hits at the transient position nor were there any strong 
gain variations (i.e. no apparent problem with the flat-fielding) at the three positions on 
the CCD. In both the sum and mode-scaled median of the three shifted images, we detect 
a faint source consistent with the centroid location (angular offset 0.04 ± 0.18 arcsec) of 
the optical transient on 1998 March 27. Lastly, we computed the point source sensitivity 
in the April 17th image by computing Gaussian-weighted photometry in 1000 random 
apertures (see discussion accompanying Table 1). Relative to this distribution, the flux 
at the location of the transient is positive and equal to 4.6-a; the probability that the 
measured flux is due to noise is 2 x 10~ 6 . All objects at the flux level of the transient are 
reliably detected in the deeper image from December, thereby providing an independent 
validation of our methodology. We conclude that indeed the transient was significantly 
detected on April 17th. We discuss the photometric calibration of the detection in Table 
1. 




27.35 Mar 1998 UT 17.25 Apr 1998 UT 18.50 Dec 1998 UT 
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Figure 2. 

The .R-band light curve of the afterglow of GRB 980326. Overlaid is a power-law 
afterglow decline summed with a bright supernova light curve at different redshifts. (Al- 
though we use as as a template the multi-band light curve of SN 1998bw (refs. 18,41), 
the bright supernova potentially associated with GRB 980425, we emphasize that the ex- 
act light-curve shapes of a supernova accompanying a GRB is not known a priori.) The 
GRB+SN model at redshift of about unity provides an adequate description of the data. 
Transient Light curve. From ref. 48 we estimate the Galactic extinction in the direction 
of the optical transient (/, b = 242°.36, 13°.04) to be E(B-V) = 0.08. Thus, assuming the 
average Galactic extinction curve (Ry = 3.1), the extinction measure is Ar = 0.22, Aj = 
0.16 mag. Plotted are the extinction corrected magnitudes (see Table 1) of the transient 
converted to the standard flux zero-point of the Cousins R filter from ref. 46) . In addition 
to our data, we include photometric detections from Groot et al. 21 and an upper-limit 
from Valdes et al. 22 (KPNO). The GRB transient flux dominates at early times, but with 
a power-law decline slope a = —2 (straight solid line). 

Supernova light curve. The supernova light curve template was constructed by spline- 
fitting the broadband spectrum measured by Galama et al. 18 of the bright supernova 
1998bw at various epochs (augmented with late-time observations of SN 1998bw by McKen- 
zie and Schaefer 41 ) and transforming back to the restframe of SN 1998bw (z = 0.0088). 
As discussed in the text, we expect the rest- frame UV emission (below 3900A) to be sup- 
pressed due to absorption by resonance lines. We assume that the UV flux declines as 
f v oc u~ s . Theoretical light curves are then constructed by red-shifting the template to 
various redshifts and determining the flux in the R (observer frame) by interpolating (or 
extrapolating, for z ^> 1). The flux normalization of the redshifted SN 1998bw curves are 
independent of the Hubble constant but are dependent upon the value of Qo and A (here 
we show the curves for fi = 0.2 and A = 0). Beyond z ph 1.3 the observed i?-band 
corresponds to restframe A <; 2900 A. As stated in the text, the spectrum in this range has 
been modeled along simple lines. Qualitatively, the peak flux derived from the SN model 
as a function of redshift and shown here agrees with the theoretical peak flux-redshift 
relation for type Type la 49 . This suggests that our adopted model for the UV spectrum is 
reasonable. 
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Figure 3 

The spectra of the transient on March 29.27 and April 23.83 1998 UT. The two spectra 
are shown at two different spectral resolutions. Starting from the top, panels 1 and 3 show 
the spectra at the two epochs at the full spectral resolution (see below for details) and 
panels 2 and 4 show the same two spectra but binned in groups of 51 channels. Panel 5 is 
the spectrum of the sky. The best-fit power law models (f^ oc v 13 ) to the binned spectra 
are shown by dashed lines; the fits were restricted the wavelength range 4500-8500 A. The 
scatter of individual channel values within each bin was used to assign relative weights to 
the median fluxes in each bin when performing the fits. 

Results. On March 29.27, we obtain (3 = -0.8 ± 0.4 and (3 = -2.8 ± 0.3 on April 23.83. 
The derived power-law indices include the correction of Galactic extinction. From the 
absence of continuum breaks in the spectrum of March 29, we can place an upper limit to 
the redshift, zot <> 2.3. 

Observing Details. Spectroscopic observations of the OT were obtained on 29 March 
1998 UT, using the Low Resolution Imaging Spectrometer (LRIS) 45 at the Keck-II 10- 
m telescope on Mauna Kea, Hawaii. We used a grating with 300 lines mm -1 blazed 
at Ablaze ~ 5000 A and a 1.0 arcsec wide slit. The effective wavelength coverage was 
A ~ 4000 — 9000 A and the instrumental resolution was ~ 12 A. Two exposures of 1800 
s each were obtained. We used Feige 34 (ref. 50) for flux calibration. The estimated 
uncertainty of the flux zero point is about 20%. Additional spectra were obtained on 
23 April 1998 UT, in photometric conditions, using the same instrument, except that the 
spectrograph slit was 1.5 arcsec wide. The effective spectral resolution for this observations 
was ~ 16 A. Three exposures of 1800 s each were obtained. For these observations we 
used HD 84937 (ref. 51) for flux calibration. The estimated zero-point uncertainty is 
about 10%. On both epochs, exposures of arc lamps were used for primary wavelength 
calibration. Night sky lines were used to correct for calibration changes due to flexure. In 
both cases, slit position angles were close to be close to the parallactic angles. Thus the 
differential slit losses were negligible. 

The spectra shown were convolved with a Gaussian with a = 5 A (i.e., less than the 
instrumental resolution) and rebinned to a common 5 A sampling. None of the apparent 
features in the spectra are real, on the basis of a careful examination of two-dimensional, 
sky-subtracted spectroscopic images: apparent emission of absorption features are all due 
to an imperfect sky subtraction noise. A sky spectrum from the April 23 observation, 
extracted in the same aperture, is shown for the comparison. These spectra are shown 
before the correction for the Galactic foreground extinction. 

Spectrophometry. In both epoch, we chose a slit position angle close to parallactic 
so that the slit would cover both the transient and a relatively bright star (R ~ 19). 
The spectroscopic i?-band magnitudes reported in table 1 were derived relative to the 
calibrated i?-band magnitude of these stars. This calibration serves to eliminate most of 
the systematics and calibration errors; that is, the spectrophotometric magnitudes were 
put on the direct CCD system, and are not based on the flux calibration of the spectra 
(which do, nevertheless, agree to twenty percent). This procedure bypasses most of the 
systematic errors in comparing our spectroscopic magnitudes with those from direct CCD 
images. 
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